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Abstract — Multiple Input Double Output (MIDO) asymmetric 
space-time codes for 4 transmit antennas and 2 receive antennas 
can be employed in the downlink from base stations to portable 
devices. Previous MIDO code constructions with low Maximum 
Lilielihood (ML) decoding complexity, full diversity and the non- 
vanishing determinant (NVD) property are mostly based on cyclic 
division algebras. 

In this paper, a new family of MIDO codes with the NVD 
property based on crossed-product algebras over Q is introduced. 
Fast decodability follows naturally from the structure of the 
codewords which consist of four generalized Alamouti blocks. 
The associated ML complexity order is the lowest known for 
full-rate MIDO codes (0{M^") instead of 0{M^^) with respect 
to the real constellation size M). Numerical simulations show 
that these codes have a performance from comparable up to IdB 
gain compared to the best known MIDO code with the same 
complexity. 

I. Introduction 

There are many wireless channel scenarios where the 
number of antennas is asymmetric, in particular, where the 
transmitter has many antennas, while the receiver, for example 
being a portable device, has few of them. MIDO channels, 
which stands for Multiple Inputs Double Outputs, refer to such 
systems. The model we will consider in this paper is a MIDO 
coherent Rayleigh fading channel, with 4 antennas at the 
transmitter and 2 antennas at the receiver, which furthermore 
has perfect channel state information at the receiver 
Since the computing power available at the receiver is typically 
very limited, the design of MIDO codes must take into account 
the decoding complexity order. A code is called fast-decodable 
if the research tree in the sphere decoding algorithm ifTOll 
can be simplified. In general, the complexity order of the 
real sphere decoding algorithm for a system with m transmit 
antennas and n receive antennas employing a space-time code 
and real signal constellations of size M is ©(M^™"). 

A. Related work 

The first low ML complexity MIDO code was proposed 
in ||2l. This code has complexity order 0{M^^) instead of 
0(M^^); it is not full-rank but still achieves good performance 
for moderate values of SNR. A full-rank MIDO code with 
0(A/^") complexity and high coding gain which is conjec- 
tured to have the non-vanishing determinant (NVD) property 
was presented in [9^1. Recently, MIDO codes with the NVD 



property based on cyclic division algebras with complexity 
order up to 0(M^°) were also constructed in iTTI l6l ITTl . 
In this paper, we consider an alternative approach and intro- 
duce a new family of 0(A/^") decoding complexity MIDO 
codes with the NVD property based on crossed-product alge- 
bras over Q. Crossed-product algebras were already used to 
construct one example of fast-decodable MIDO code in |7|; 
however, this example is based on puncturing a 4 x 4 full-rate 
space-time code and its performance is not very good. The 
constructions presented here are tailored for the 4x2 case 
and do not require any puncturing. 

B. Alamouti-like structures 

Let 0* denote the complex conjugation for a scalar, and 
the Hermitian transpose for a vector or a matrix. Recall that 
an Alamouti block code is given by 

( 2/1 -y*2\ 
\V2 yl ) ' 

with 2/1,2/2 G C. It has the property that its columns are 
orthonormal. This property allows fast decoding, and con- 
sequently many of the attempts to construct fast decodable 
space-time codes have tried to mimic it. In this paper we will 
consider generalized Alamouti codewords of the form 

{ I' ,aeC (1) 

\a*y2 yt J 

which will make the columns orthogonal (furthermore or- 
thonormal if |a| = 1). Note that this definition is different 
from the one in (T\. Note also that if we have a matrix of the 
form 

-J*A,aeR, a>0, (2) 

then by multiplying the second column by l/y/a, and the 
second row by y/a, we get, without changing the matrix 
determinant, 

( P , (3) 

Vv"2/2 yl J 

a particular case of ([T]i when a is real. As far as decoding 
is concerned, it is enough to ask for the columns to be 
orthogonal; the orthonormality does not improve the decoding 



complexity, but rather the performance by ensuring that the 
energy is balanced across time and antennas. 

Our goal in this paper is to construct a code carrying 8 
complex information symbols of the form 

'A C 
B D 



1) 



(4) 



where the 2x2 blocks A and D are generalized Alamouti 
codes of the form ([T]) (it will be discussed in Section IV why 
the focus is on the blocks A and D), preferably with columns 
as close to orthonormal as possible. This means that the energy 
of the symbols might not be balanced, as in ([3]l. 

II. The Framework of Crossed Product Algebras 

The incentive to consider biquadratic crossed product al- 
gebras as underlying algebraic structure to construct fast 
decodable MIDO codes is two-fold: first, we will see below 
that the representation of these algebras naturally gives rise 
to codewords of the form Q, and furthermore, as is the case 
with traditional space-time coding using division algebras, a 
codebook with full diversity is obtained from division crossed 
product algebras. 

A. Crossed product algebras of degree 4 

We will consider as in 11] a crossed product algebra A = 
{L/K, a, b, u) over the biquadratic extension L/K, where L = 
K{Vd, Vd'), and 

(a) = Gal{K{VI)/K), (r) = Gal{K{Vd)/K). 



Such an algebra is of the form A = L 

o2 



eL®fL® efL, where 
a e K{Vd), P = b e K{Vd'), xe = ea{x) Wx G L, 
xf = fT{x) Vx G L, and fe — efu for u a non-zero element 
of L such that ua{u) — a/T{a), ut{u) = a{b)/b. Elements 
of A admit the following matrix representation: 



Xi 
V ^3 

with xo,xi 



aa{xi) 
a{xo) 
T{a)ua{x3) 
ua{x2) 

X2,X3 G L. 



6t(x2) abT{u)aT(x3) \ 

6t(x3) bT{u)(TT{x2) 

T(xa) T{a)(TT{xi) 

t(xi) crT{xo) ) 



(5) 



As shown in HI, u is such that Ni^/x{u) — 1, and suitable a 
and b are determined by the choice of u in the following way: 



'k^/d 

l{l + ua{u)) 



if ua{u) - 
otherwise, 



-1, keK; 
I G K. 



(6) 



In the latter case, in order to have a non-degenerate algebra, 
ua{u) should belong to K{y/d) \ K. Similarly, 



if ut{u) = 
otherwise, 



-1; k' e K 
I' G K, 



(7) 



l+'iir (ti) 

with ut{u) G K{^/d/) \ K. However, checking whether the 
resulting crossed product algebra ^ is a division algebra only 
depends on the choice of u: 

Theorem 1. STj Let K be a number field, and let A = 
{L/K,a,b,u) be a crossed product algebra. Then A is a 
division algebra if and only if 



K 



is a division algebra and ua{u) = — 1, 



2) ^ — — — K{s/d}/K( — '^^J ^ division algebra and 
ua{u) 7^ —1. 

There is a similar equivalent formulation depending on 
whether ut{u) = — 1. Since we need 8 complex symbols, 
it is enough to consider a crossed product algebra A over a 
biquadratic extension L of K ~ Q^. Such an algebra is of index 
4, thus we can encode 16 real symbols, that is 8 complex ones. 
In order to obtain a matrix of the form Q from (|5]l. 



A 



Xq aa{xi) 
xi a{xo) 



D 



't{xo) T{a)aT{xi) 
t{xi) (tt{xq) 



should be generalized Alamouti blocks. Alternatively, by 
swapping the second and third rows and the second and third 
columns in (j5]l, we get 



A 



Xo bT{x2) 

X2 r{xo)l ' 



D 



a{xo) bT{u)aT{x2) 
ua(x2) ctt^xq) 



Finally, by swapping the second and fourth rows and 
columns, we obtain 



A = 



xq abT{u)(jT{x3 

2^3 0-^(2:0) 



D = 



t{xo) T{a)ua{x3)\ 
brixs) a{xo) J 
(8) 

Q{Vd) and 



There are three possible choices for K{y/d) 

• both are imaginary quadratic fields, 

• both are real quadratic fields, 

• one is a real quadratic field, the other is an imaginary 
field. 

Since typical signal constellations such as QAM are encoded 
using Q{i), we will assume that one of the two quadratic fields 
is and will thus not consider the case of two real quadratic 
fields. Due to the lack of space, we will focus on the case of 
two imaginary quadratic fields, where err acts as the complex 
conjugation, and choose codewords of the form ([8]l. 

B. Crossed product algebras over imaginary fields 

Consider now the case where d = —c < 0,d' = —c' < 0, 
that is L = Q(-\/--c, -\/~c'). We are mostly interested in the 
case when c' = 1, but nevertheless will show later that the 
construction is also possible for other values of c'. 



The Galois group of Q{i 
by (a), resp. (t) with a{i) = 
element of L is of the form x 



resp. (1, 

—i, t(^— c) = 
: ai+a2i + a3 



is denoted 
c. Every 



Oi J 0,2 J 0-31 0,4 G Q, and we extend a and r to i so as to get 

(t{x) = fli — + a^\/—c — 041%/— c, 
T(a::) = Oi + — a3-\/^ — a^i\f^^^ 
(jr{x) — ax — a-ii — a3\/— c + 042^/^ = x*. 

We need to obtain two properties: 

• the Alamouti-like block structure for fast-decodability, 

• the algebra should preferably be a division algebra to 
guarantee a good behaviour at high SNR. 

Let us first exploit the property aT(x) — x* to construct 
a code with an Alamouti block structure by swapping the 
second and fourth row and the second and fourth column of 
the representation (j5]l, as already mentioned in ([8]): 

/ xq abT{u)aT{x^) hT{x2) aa{xi) \ 
xs, (7t{xo) t{xi) U(t{x2) 

X2 T{a)(jT{xi) t{xq) T{a)u(T{xz) 

\ Xi bT{u)(TT{x2) bT{x3) Cr(xo) / 

Since complex conjugation commutes with cr and t, we can 
rewrite it as 

/ Xq abT{u)x'!^ hT{x2) aT{xi)* \ 

X3 Xl t{xi) Ut{x2)* 

X2 T{a)xl t{xo) T{a)uT{x3)* 
\ xi bT{u)x^ bT{x3) t(xo)* / 

We are now left with the choice of a, b and u. From (|2]), we 
have as condition on these parameters that 

abT{u) e M, abT{u) < 0. (10) 

To obtain a suitable crossed product algebra A, we need to 
choose the element u (with Ni^/}^{u) = 1), and take a,b such 
that 




and 



c, k eQ 

ua{u)), I e 



k'i, k' eQ 



if uij(u) = —1; 
otherwise, ua{u) ^ K, 



if ut{u) = 
otherwise, 



-1; 

ut(u) ^ K. 



(11) 



(12) 



l+iir(u) 

For such an algebra to be a division algebra, when ua{u) = 
— 1, we have by Theorem [T] that it is enough to check whether 
^^^^ is. This in turn is equivalent to see whether c — 
-/VQ(i)/Q(s) for some s e Q{i). Since iVQ(i)/Q(s) = sj + si 
for s — si +is2, si,S2 G Q, we finally need to check whether 
c can be written as a sum of two squares. We also consider 
the case where c' ~ 2, which is the next smallest c' after 
c' = 1. Similarly in this case when ua{u) — —1, we need to 
check whether (^^15^) is ^ division algebra, that is, whether 
c — s\ + 2s\, si, S2 G Q- Recall that c can be written as the 
sum of two rational squares if and only if all its odd prime 
factors which are congruent to 3 (mod 4) occur to an even 
exponent; similarly, c can be written in the form s\ + 2s\ if 
and only if its odd prime factors which are congruent to 5 
or 7 (mod 8) occur to an even exponent. The first smallest 



possible values for c are listed below, for both c' 
c' = 2 0. 



1 and 



(c,-c') 


division alg 


ebra 


(c,-c') 


division alj 


ebra 


(2,-1) 


no (2 = 1 


f 1) 


(2,-2) 


no (2 = 


+ 2) 


(3,-1) 


yes 




(3,-2) 


no (3 = 1 


+ 2) 


(5,-1) 


no (5 = 1 


f 4) 


(5,-2) 


yes 




(6,-1) 


yes 




(6,-2) 


no (6 4 


+ 2) 


(7,-1) 


yes 




(7,-2) 


yes 




(10,-1) 


no (10 = 9 


+ 1) 


(10,-2) 


yes 




(11,-1) 


yes 




(11,-2) 


no (11 = 9 


+ 2) 


(13,-1) 


no (13 = 9 


+ 4) 


(13,-2) 


yes 





The condition when u(j{u) 7^ —1 is a priori less systematic 
to check, since it relies on the trace of ua{u), though we can 
look for an element u such that the trace of ua{u) is c — 2, 
where c is such that (c, — c') is a division algebra. 
We finally discuss briefly how to find u with Ni^/j^{u) = 1. A 
natural choice is to start by taking u a unit in L. Recall that by 
Dirichlet's unit theorem, the units of an algebraic number field 
i are a multiplicative group generated by a set of fundamental 
units. The number of fundamental units is ri+r2 — 1, where ri 
is the number of real embeddings of L, and r2 is the number 
of pairs of complex embeddings. 

In the case L — Q{y/^, \/—c'), we have ri = and r2 — 2, 
therefore there is only one fundamental unit. 

III. Code Constructions 

A. A generic construction 

Suppose that there exists u in L with Ni^/x{u) = 1, and 
that the corresponding a and b as defined by (111 and (12i 
satisfy the following conditions: 

U<7{u) ^ —I, MT(?i) = e G {— 1, i, — «}, 
abT{u) G M, abT{u) < 0. 

We now set 

a = —abT{u) > 0. 

Observe that the first condition implies that a = iy/ck, and 
take k = I. Since b G K{y/^) = Q(i), a{b) = b* and the 
second condition, with I' = 1, implies that e = ut{u) 
We then have 



(13) 



b ■ 



bT(u) 



a 
a 



u 



t{u) 
T{a)u — - 



eab 
a 



ebi-Jc i^/cb* 



a 



a 



/ cu ■ 



cbe cb* 
a a ' 

where the second equality uses the expression computed for 
t{u) above. Replacing in the expression (j9|, we find that a 
codeword is of the form 



/ 2^0 




bT{x2) 


i^/cT{xi)* \ 




Xq 


t{xi) 


--^r{x2y 


X2 


—i\fcx\ 


t{xo) 


-'^rix.r 


V Xl 


ia * 
v^^2 


brixs) 


rixoT J 



Dividing the first row by -y/a and multiplying the first column 
by -y/a, and further multiplying the fourth column by and 



dividinj 


J the fourth row 


by yields: 




/ 








V 




t{xx) 




V 


—i\fcx\ 


t{xq) -- 






fcxx i^/cix\ 





irixiY \ 
r{xo)* ) 



(14) 



We have thus obtained a codeword composed of four gener- 
alized Alamouti blocks; 



/ ^1 



Zb 

z& 

Zl 
Z8 



6 



—iz 



8 



We now provide examples of such code constructions, with 
values of c which give a division algebra, namely c' = 1 and 
c = 3,6 and 11. In the following, the fundamental units have 
been computed using the KASH software |4|. 

Example 1 (L = Q{i,y/3)). Let c' = 1, c = 3. Then 
L = Q{i, a/3) = Q(Ci2)- The fundamental unit of L is u = 
(i±^) (V^-1). We will choose u = t{v) = {-V3-1). 
We have u<j{u) = — 1, ut{u) = —i. Then a = \/— 3 and 
& = We have a6r(u) = \/3(i^) < 0, so the 

conditions ^T3\ are satisfied. 

Example 2 (i = Q(j, a/6)). Let c' = 1, c = 6. We can choose 
as u the fundamental unit u — (l + «) ^ Ij- We have 



ua{u) — —1, ut(w) = — i. Therefore a 
in the previous case, abT{u) — 



Example 3 {L = Q{i, 
choose as u the unit u = 



\/2 



As 



11)). Let c' = 1, c = 11. We can 

(-3-\/n). We have ua{u) = 



-1, ut{u) = —i. Therefore a = i\/Ti, b = Again, we 



have abT{u) = vH 



3-\/ll 



< 0. 



We conclude by giving an example with c' = 2 instead of 
c' = 1, with c = 5 to get a division algebra. 

Example 4 {L = Q{iy/2,iy/5)). Let c' = 2, c = 5. The 
fundamental unit is u = 3 — a/IO- We have ua{u) = —1, 
ut{u) = —1, a = ia/S, b = za/2, and abT{u) = — a/T0(3 + 

Vio) < 0. 

Remark 1. Ideally, the parameters u, a and h ought to be 
of complex norm 1 in order to have good energy efficiency 
|[r|. Unfortunately, this does not seem to be possible. Indeed, 
in the case of two imaginary quadratic subfields K{y/—c) and 
K{\/~c'), if a and b have complex norm 1, since the automor- 
phisms a and r act as the complex conjugation respectively 
on K{^/^^ and K{\/—cJ), we have aria) — — 1 and 
b(j{h) = \b\ = 1. As a consequence, the condition in Theorem 
[Tjcannot hold and the crossed-product algebra is not a division 
algebra, since Vq G Q, the quaternion algebra ' — 



a division algebra. Clearly, the case of real quadratic subfields 
is also hopeless because they do not contain any non-trivial 
roots of unity. 

We finally give an example of a crossed product algebra 
which is not a division algebra, but provides good shaping: 

Example 5 (i = Q{i,V3)). Let c' = 1, c = 3 as in 
Example [1] and choose u = ( = (12- We have uaiu) = C^, 
ut{u) = -1. Then a 1 + C** and 6 = i. We have 

abT{u) = —1, so the conditions (13i are satisfied. Since u, a 



and b have complex norm 1, this code provides a very good 
shaping, although unfortunately it is not full-diversity. 

B. Code optimization 

As we have seen, in order to obtain an Alamouti-like 
structure, the price to pay is to unbalance the energy in the 
codewords. We now consider the code C in Example [T] and 
discuss some changes to improve its performance. 

Choice of the basis of L: Recall that an element of L 
is of the form x — ai + a2\J —d + asa/^ + a/^\J —d \f^, 
01,02,03,04 e Q. Thus Examples 1-3 and 5 allow QAM 
encoding, since: 

X = ai+a2i+a3A/^+a3A/^7; = (ai+a2z) + A/^(a3+ia4) 

where both a\ + 02* and 03 + 104 are in Q(i). 
In particular, to encode QAM symbols, we can choose any 
basis of L as a vector space over Q(i). We will consider two 
bases, which may not coincide: the basis B\ of the ring of 
integers Ol over Z(i), and the basis B2 — {1,\/^}. B\ 
corresponds to a denser lattice, but we will see that B2 is 
more convenient for fast decodability, since it is composed by 
one purely real and one purely imaginary element. 
Let C2 be the version of the code C employing the basis 
B2 = {l,A/3i}. Its drawback is that the basis vectors have 
unbalanced norms. In order to improve the shaping, we also 
consider the code C3 with basis S3 = {2, a/Si}. 

Renormalization of the parameters a and b: Since a and 
b are only defined up to a rational constant in ^ and (j7|, we 
can choose the normalization in such a way that their complex 
norm is close to 1, for example we can choose k — I' — ^ in 
Example [T] (code C4). However, this renormalization affects 
the minimum determinant of the code; indeed, if k and I are 
not integers, the code is not contained in the natural order 
of the crossed-product algebra. A further improvement can be 
obtained by multiplying the block C and dividing the block 
B by in J4J1 by |a|^; this does not affect the determinant. 

IV. Fast Decodability 

In this section we briefly review some facts about decoding 
complexity. Consider a linear dispersion code encoding K 
real symbols ai, . . . , qk in a constellation S such that each 
codeword is of the form X — J2iLi Mo-i with generator 
matrices Ai, . . . , A^- The code is called conditionally g-group 
decodable |8| if there exists a partition of {1, ... , into g+\ 
disjoint subsets Fi, . 

bij = \\AiA^ + Ar, 



Fg, F'^ such that 



I \\f 



e Fi,VTO e rj,i 7^ j. 



4-QAM 




10 



-B- Code A4 (complexity 10) 
-*- CodeC4 
-I- CodeC2 
-e- Codecs 

Code C5 (not full diversity) | 



10 



12 14 16 

SNR (dB) 



20 



Figure 1. Performance comparison of the proposed codes using 4-QAM 
constellations with the version of the A4 code with A/^" decoding complexity. 



In this case, the sphere decoding complexity order reduces to 
^^|r^|maxi<.<,|r,|^ where M is the size of S. 
For codes of the form ([T4]l, using the basis B2, one can check 
that the Hurwitz -Radon matrix B = [hij) L3j has the shape 



where t denotes any possibly nonzero symbol. Clearly, this 
code is conditionally 4-group decodable and has complexity 
order 10. In fact, in order to decode one can Ust all the possible 
values for the variables {xg, a;io, . . . , xiq] and then minimize 
separately the Euclidean distance over the pairs of variables 
{xi,X2}, {0:3,2:4}, {2:5,2:6} and {xt,x^}. 
When using the basis Bi, one can see that the complexity 
order is 12 (details are omitted for lack of space). 

V. Simulations 

Figure [T] shows the performance of the proposed codes 
using 4-QAM, compared to the 0(M^°) decoding complexity 
version of the "yl4 code" in 1 12l at the same spectral efficiency. 
The code C2 loses about 1.6 dB with respect to the A^^ code 
at the WER of 10^"'. Thanks to its more balanced basis, 
the code C3 improves the performance by 1 dB. Finally, the 
renormalized version C4 has better performance than the ^4 
code in the low SNR regime (up to SNR=18) but does not work 
so well at hish SNR due to its small minimum determinant. 



Due to its excellent shaping, the code C5 from Example 
|5] performs surprisingly well at low SNR, and in fact it 
outperforms the code by 0.9 dB at the WER of 10^*, even 
though the error rate will eventually be worse in the high SNR 
regime since the code is not full rank. 

VI. Conclusions 

We proposed a new family of full-rate, non-vanishing de- 
terminant 4x4 MIDO codes based on cross-product algebras 
over Q with ML decoding complexity order 0{M^^). 
Simulation results show a performance from similar up to a 1 
dB gain compared to the best previously known code with the 
same complexity order 1 12 1. While the latter requires real PAM 
signal constellations, our codes further have the advantage of 
being suitable for QAM complex modulation. 
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